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The  Departments  of  Aerospace  Engineering  and  Mechanical  Engineering  at  the  University 
of  Southern  California  acknowledge  support  from  the  Office  of  Naval  Research  under  contract 
number  NOOOl 4-82- K. -0084  for  the  period  December  1,  1981  through  September  30,  1988.  The 
overall  objective  of  the  research  was  (and  continues  to  be)  to  understand  through  laboratory 
experiments,  theoretical  analyses,  and  comparisons  with  field  data,  the  processes  by  which  turbulence 
in  the  world’s  oceans  is  generated,  evolves  and  decays.  Laboratory  studies  are  judiciously  chosen  to 
illustrate  fundamental  processes.  Efforts  are  made  to  relate  the  results  of  laboratory  studies  to  the 
ocean  environment. 


During  the  seven  years  of  the  contract  period,  a  total  of  56  published  articles  have  appeared 
in  journals  or  conference  proceedings.  In  addition,  two  reports  have  been  prepared  and  given  wide 
circulation  within  the  research  community.  Highiights  of  work  during  the  contract  period  are  the 
following. 

I.  1.  Comprehensive  studies  of  wave  motions  in  rotating  flows,  stratified  flows,  and 
_ rotating /stratified  npws  ;^ 

Wave  motions  are  important  in  transporting  momentum  and  energy  within  the  oceans  and  atmosphere. 
The  laboratory  studies  have  served  to  deepen  and  extend  the  understanding  of  such  waves.  We  were/ 
the  first  to  describe  solitary  wave  motions  on  concentrated  vortex  cores.  A  series  of  fundamental 
experiments  have  embraced  a  wide  range  of  applications  including:  the  dynamics  of  severe  storms 
in  the  atmosphere,  and  the  quantized  vortex  motion  in  liquid  Helium.  Solitary  Kelvin  waves  —  long 
wave  motions  present  in  ail  the  world’s  oceans  --  were  first  studied  in  the  laboratory  here.  The 
mechanism  for  tidal  generation  of  internal  gravity  waves  (common  features  of  all  continental  shelf 
regions)  was  also  first  proposed  and  studied  here. 

II.  SPvnamics  of  gravity  currents  . 

The  first^udies  of  axisymmetric  internal  gravity  currents  were  carried,  out  at  USC,  as  well  as  the 
first  studie’s  of  gravity  currents  having  a  variable  rate  of  supply.  Both  th^se  circumstances  are  more 
typical  of  naturally  occurring  gravity  currents,  and  are  therefore -important. 


III.  ■  Effect  of  stratification  upon  the  structure  of  oceanic  turbulence 


Turbulence  in  the  ocean  must  is  a  result  of  the  stratification  of  the  medium  ^Continuously  lift 
heavy  fluidtfrom  below)Jand  depress  lighter  flui^tfrom  above).*  A  consequence  is  that  vertical  scales 
of  turbulence  are  limited  in  extent.  We  have  demonstrated  this  limiting  scale  in  several  dramatic 
experiments. 


The  turbulent  eddies  continue  to  evolve  and  spread  in  the  horizontal  direction,  producing  a  collection 
of  pancake-shaped  structures.  This  two  dimensional  (or  quasi-twp  dimensional)  turbulent  eddy  field 
is  characteristic  of  the  ocean  at  scales  from  hundreds  of  meters  to  hundreds  of  kilometers.  Much  of 
our  recent  effort  has  gone  into  undentanding  the  evolution  of  this  complex  field. 

IV..f  V  The  common  features  of  frontal  motions /coastal  upwelling 

The  major  fisheries  of  the  world  occur  in  regions  of  upwelling.  The  two  ingredients  are  a  coast  line 
adjacent  to  deep  water  and  an  along-shore  wind  field.  Ekman  flux  in  the  surface  boundary  layer  -? 


produces  upwelling  aJong  the  coast.  Studies  in  our  laboratory  show  the  propagation  of  waves  along 
the  upwelling  front,  and  often  subsequent  instability  and  break-up  into  off-shore  Jets.  The 
resulting  structure  is  highly  dependent  upon  the  local  topography.  These  frontal  instabilities  are 
one  source  of  the  two  dimensional  turbulent  eddies  mentioned  in  III.  The  laboratory  studies  yield 
scaling  laws  which  are  applicable  to  the  oceanic  counterparts.  Fractal  geometry  is  employed  here  for 
the  first  time,  as  a  means  of  clarifying  the  motion  and  to  characterize  the  degree  of  mixing. 

V.  Fundamental  structure  of  turbulent  flows  ^ 

In  addition  to  studies  of  oceanic  turbulence,  we  also  investigate  fundamental  issues  on  the  nature  of 
turbulence.  These  experiments  are  more  easily  accomplished  in  homogeneous  (unstratified)  Hows, 
although  the  results  have  broad  application  to  all  turbulent  flows.  We  were  the  first  to  emphasize  the 
geometrical  characterization  of  vortex  motions  in  free  shear  flows.  We  have  identified  certain  defect 
features  which  seem  to  be  fundamental  to  ail  unbounded  shear  flows,  and  have  devised  a  way  co 
produce  and  study  these  defects  artificially. 

VI.  .^  Dev.glODment  and  implementation  Qf_ image  processing  algorithms  ; 

The  two  published  reports  deal  with  the  general  problem  of^article  tracking  in  fluids.  The 
<*panicles^ay  be  floats  in  the  ocean  or  neutrally  bouyant  particles  added,  to  a  laboratory  flow  for 
the  purpose  of  visualizing  the  fluid  motionf^  Several  tracking  and  grid  interpolation  schemes  were 
evaluated  in  the  report  and  form  the  basis  of  a  free  software  package  which  is  included  with  the 
report  on  a  inch  floppy  disk. 


/ 
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Turbulence  and  Waves  in  a  Rotating,  Turbulent  Flow",  with  F.  Browand  and  E.J.  Hopfinger. 

J.  Fluid  MedL.  121.  PP-  505-534,  1982. 

"Vortex  Solitary  Waves  in  a  Rotating  Turbulent  Flow",  with  F.  Browand  and  E.J.  Hopfinger. 
Nature.  295.  February  1982. 

The  Dynamics  of  Confined  Vortices",  with  T.  Maxworthy,  M.  Escudier  and  J.  Bornstein. 
Proc.  Rov.  ^c..  Ser.  .A.  p.  335,  1982. 

"Observations  of  Inertial  Waves  in  a  Homogeneous.  Rotating  Fluid",  with  T.  Maxworthy  and 

K. .  Heikes,  J.  Ruid  Mech..  121,  P.  319,  1982. 

"Experiments  on  Solitary,  Internal  Kelvin  Waves".  T.  Maxworthy,  J.  Fluid  Mech..  ii2.  PP- 
365-383,  1983. 

"On  the  Generation  and  Evolution  of  Internal  Gravity  Waves",  with  T.  Maxworthy  and  F. 
Lansing,  J.  Fluid  Mech..  p.  124,  1984. 

"Wave  Motions  on  Vortex  Cores",  with  T.  Maxworthy,  E.J.  Hopfinger  and  L.G.  Redekopp, 
J.  Huid  Mech..  lii.  pp.  141-165,  1985. 

"Solitary  Waves  on  Density  Interfaces’,  T.  Maxworthy,  Proc.  of  Conference  on  "Waves  on 
Interfaces".  Madison,  Wisconsin,  October  1982. 

"Dynamics  of  the  Great  Red  Spot  Inferred  from  Voyager  I  and  II  Images",  with  T.  Maxworthy 
and  J.L.  Mitchell,  submined  to  ICARUS. 

"Waves  on  Vortex  Cores  and  Their  Relationship  to  Vortex  Breakdown",  with  T.  Maxworthy, 
M.  Mory  and  E.J.  Hopfinger,  Proc.  .aGaRD  Mta.  on  Vortex  Flows.  Brussels,  1983. 

"Internal  Wave  Generation  and  Propagation  in  a  Rotating-Stratified  Fluid",  with  T. 
Maxworthy,  G.  Chabert  d’Hieres  and  H.  Didelle,  J.  Geoohvs.  Res..  89.  pp.  6383-  SB6,  SW 

The  Generation  of  Kelvin  Waves  by  Tidal  Flow  over  Bottom  Topography  in  a  Narrow 
Channel",  with  T.  Maxworthy  and  Y.N.  Liu,  J.  Geophvs.  Res..  2fl.  PP-  7235-7242,  1985, 

"Wave  Motions  on  Concentrated,  Columnar  Vortices",  T.  .Maxworthy,  M.  Kruskal  Festschrift. 
Special  issue  of  Phvsica  D.  ij,  p.  149,  1986. 

"Waves  on  Vortex  Cores",  T.  Maxworthy,  Proc.  lUTAM  Svmo.  on  "Fundamental  Aspects  of 
Vortex  Motion".  Tokyo,  August  31-September  4,  1987,  (Invited  lecture). 


Dynamics  of  gravity  currents 

The  Viscous  Spreading  of  Plane  and  Axisymmetric  Gravity  Currents",  with  T.  Maxworthy 
and  N.  Didden,  J.  Fluid  Mech..  liL,  P-  27,  1982. 


The  Dynamics  of  Double-Diffusive  Gravity  Currents",  T,  Maxworihy, 
pp.  259-282,  1983. 
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“Gravity  Currents  with  Variable  Inflow",  T.  Maxworthy.  J.  Fluid  Mech..  128.  pp.  247-257. 
1983. 


III.  Effect  of  stratification  uoon  the  structure  of  oceanic  turbulence 

1.  "Selective  Withdrawal  and  Spin-up  of  a  Rotating.  Stratified  Fluid",  with  T.  .Maxworthy  and 

S.  Monismith.  accepted  for  publication  in  J.  Fluid  Mech. 

2.  "Fractal  Geometry  in  the  Energy  and  Vorticity  Fields  in  Decaying  2D  Grid  Turbulence’,  with 

T.  .Maxworthy.  G.  Spedding  and  Ph.  Caperan.  submitted  to  J.  Geophvs.  Res. 

3.  “Interaction  de  Tourbillons  Bidimensionelle  de  Meme  Signe:  Etude  Experimentale  et 
Numerique",  with  T.  Maxworthy,  Ph.  Caperan.  J.  Verron  and  E.J.  Hopfinger,  Proc.  of 
ColloQuium  "Dvnamiaue  de  fluides  Geophvsioues  et  AstrophvsiQues".  Grenoble,  France. 
September,  1987. 

4.  The  Stratified  Mixing  Layer",  with  F.  Browand,  G.A.  Lawrence  and  J.C.  Lasheras,  6th 
Conference  on  Turbulent  Shear  Flows,  Toulouse,  France,  August  1987. 

5.  "The  Evolution  of  Grid -Generated  Turbulence  in  the  Neighborhood  of  a  Density  Step",  with 

F.  Browand,  R.  Amen  and  T.  Maxworthy,  for  the  lAHR  3rd  Symposium  on  Stratified  Flow, 
Cal.  Tech.,  Pasadena,  CA,  Feb.  3-5,  1987,  published  Proceedings. 

6.  "Shear  Instabilities  in  Stratified  Flow",  with  F.  Browand,  G.A.  Lawrence  and  J.C.  Lasheras, 
for  the  lAHR  3rd  Symposium  on  Stratified  Flow,  Cal.  Tech.,  Pasadena,  CA,  Feb.  3-5,  1987, 
published  Proceedings. 

7.  "The  Behavior  of  a  Turbulent  Front  in  a  Stratified  Fluid:  Experiments  with  an  Oscillating 
Grid",  with  F.  Browand,  D.  Guyomarand  S.-C.  Yoon,  for  the  lUTAM  Symposium  on  Mixing 

.  in  Stratified  Fluids,  Margaret  River,  Western  Australia,  August  25-29,  1985,  Also  discussed 
at  APS,  Division  of  Fluid  Dynamics  Meeting,  Tucson.  .Arizona,  Nov.  1985.  J.  Geophvs.  Res.. 
Vol.^,  pp.  5329-5341,  1987. 

8.  "Quasi  2-D  Turbulence  and  Vortices  in  a  Stratified  Fluid",  with  T.  Maxworthy,  Ph.  Caperan 
and  G.R.  Spedding,  3rd  Int.  Conf.  on  Stratified  Flows.  Caltech  1987. 

9.  "Collapse  of  a  Turbulent  Front  in  a  Stratified  Fluid.  Part  1.  Nominally  Two-Dimensional 
Evolution  in  a  Narrow  Tank",  with  T.  Maxworthy,  Y.N.  Liu  and  G.R.  Spedding,  J.  Geophvs. 
Res..  22,  pp.  5427-5433,  1987. 

10.  "The  Structure  of  Meso-Scale  Baroclinic  Turbulence",  with  T,  Maxworthy,  S.  Narimousa  and 

G. R.  Spedding,  Sixth  Conference  on  Turbulent  Shear  Flows.  Toulouse.  France,  August  1987. 

11.  "An  Experimental  Investigation  of  the  Coalescence  of  Two-Dimensional,  Finite-Core 
Vortices",  with  T.  Maxworthy  and  Ph.  Caperan,  submitted  to  J.  P.  O. 

12.  "Experiments  in  Two-Dimensional  Turbulence;  the  Interaction  and  Merging  of  Vortices", 
with  T.  Maxworthy,  Ph.  Caperan  and  G.R.  Spedding,  Proc.  Inter.  Conf.  on  Fluid  Mechanics. 
Beijing,  China,  July  1987. 

13.  "The  Effects  of  Rotation  on  Selective  Withdrawal  from  a  Density-Stratified  Reservoir",  with 
T.  Maxworthy  and  S.  Monismith,  3rd  Int.  Svmp.  on  Stratified  Flows.  Caltech  1987. 

"A  Note  on  Turbulent  Mixing  Across  a  Density  Interface  in  the  Presence  of  Rotation",  T 
Maxworthy,  J.  Phvs.  Ocean.,  ifi,  p.  1136,  1986. 


14. 


15. 
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"Natural  Flow  Visualisation  in  Lake  Eyre",  with  T.  .Maxworthy,  J.  Bye.  R.  Byron-Scott  and 
R.  Nunes,  EOS  67.  No.  5.  South  Australia.  Feb.  4,  1986. 

16.  The  Fluid -Dynamics  of  Natural  Phenomena".  T.  Maxworthy,  lOth  U.  S.  National  Congress 
of  ADolied  Mechanics,  to  appear  in  A.S.M.E.  Trans..  1986.  Invited  lecture. 

17.  "Differential  .Mixing  in  a  Stratified  Fluid*,  with  T.  .Maxworthy,  and  S.  Monismith,  J.  Fluid 
Mech..  189.  pp.  571-598,  1988. 

18.  "Differential  Deepening  of  a  Mixed  Layer  in  a  Stratified  Fluid",  T.  Maxworthy,  I.A.H.R. 
2Jst  Congress,  .Melbourne,  Australia,  August  1985. 

19.  The  Inhibition  of  Vertical  Turbulent  Scale"  with  F.  Browand  and  E.J.  Hopfinger,  IMIA 
Conference  Proceedings  on  Models  of  Turbulence  and  Diffusion  in  Stably  Stratified  Regions 
of  the  Natural  Environment,  Cambridge,  England,  March  14-15,  1983. 


IV.  The  common  features  of  frontal  motions /coastal  upwelling 

1.  "Laboratory  Models  of  Oceanic  Fronts",  T.  Maxworthy,  4th  Asian  Conference  on  Fluid 
Mechanics,  Hong  K.ong,  September  1989.  Also  submitted  to  J.  of  Marine  Research. 

2.  "Fingering  Instability  at  High  Flow  Rates  in  a  Hele-Shaw  Cell",  with  T.  Maxworthy,  G.R. 
Spedding  and  J.  Wiggert,  in  preparation  for  J.  Fluid  Mech. 

3.  "An  Experimental  Study  of  Interface  Instability  in  a  Hele-Shaw  Cell",  T.  Maxworthy,  1st 
Caribbean  Conf.  on  fluid  Dynamics,  Univ.  of  the  West  Indies,  Trinidad,  January  1989.  Also 
submitted  to  Phvs.  Rev.  A 

4.  The  Non-linear  Growth  of  a  Gravitationally  Unstable  Interface  in  a  Hele-Shaw  Cell",  T. 
Maxworthy,  J.  Fluid  Mech..  177,  pp.  207-232,  1987.  Also  a  note  in  Phvs.  Fluids.  p. 
2637,  1985  (1st  Prize  A.P.S.-D.F.D.  Photo  Competition). 

5.  "Coastal  Upwelling  on  a  Sloping  Bottom;  The  formation  of  jets,  plumes  and  pinched-off 
cyclones",  with  T.  Maxworthy  and  S.  Narimousa,  J,  Fluid  Mech..  JT^,  pp.  169-190,  1987. 

6.  "Coastal  Upwelling  Driven  by  a  Discontinuous  Stress  Distribution  in  the  Presence  of  Bottom 
Topography",  with  T.  Maxworthy  and  S.  Narimousa,  J.P.O..  16.  PP.  2072-2083,  1986. 

7.  The  Laboratory  Modelling  of  Coastal  Upwelling",  with  T.  Maxworthy  and  S.  Narimousa, 
Oceanography  Report  Section  of  EOS  67.  No.  1,  January  7,  1986. 

8.  "Bubble  Formation,  Motion  and  Interaction  in  a  Hele-Shaw  Cell",  T.  .Maxworthy,  J.  Fluid 
Mech.  G.l.  Tavlor  Special  Volume.  December  1986. 

9.  "Experiments  in  a  Hele-Shaw  Cell",  T.  Maxworthy,  in  Experiments  in  Fluid  Dynamics,  ed. 
R.C.  Granger,  1986. 

10.  "Experiments  and  Observations  on  the  Effects  of  Bottom  Topography  on  Coastal  Upwelling, 
with  T.  Maxworthy  and  S.  Narimousa,  Proc.  of  Conf.  on  Vertical  Motion  in  the  Tropical 
Oceans.  Paris,  May  1985. 

11.  "The  Formation  of  Tidal-Mixing  Fronts",  T.  Maxworthy,  Ocean  Modelling.  No.  55,  March 
1984. 
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12.  *A  Two-Layer  Model  of  Shear  Driven  Coastal  Upwelling  in  the  Presence  of  Bottom 
Topography",  with  T.  Maxworthy  and  S.  Narimousa,  J.  Fluid  Mech..  159.  pp.  503-531,  1985. 


V.  Fundamental  structure  of  turbulent  Hows 

1.  "Pattern  Evolution  in  the  Two-Dimensional  Mixing  Layer",  with  F.  Browand  and  S.  Prost- 
E)omasky,  NATO  Workshop,  New  Trends  in  Non-linear  Dynamics  and  Pattern  Forming 
Phenomena;  The  Geometry  of  Non-equilibrium,  Cargese,  France,  August  2-12,  1988, 
Proceedings  to  be  publis|ed  in  Phvsica  D. 

2.  ".A  Technique  for  Acoustic  Excitation  of  Separated  Shear  Flows:  Preliminary  Results",  with 
F.  Browand  and  S.  Prost-Domasky,  Proceedings  of  the  .aSME  Winter  Meeting.  Chicago, 
Illinois,  November  28-Dccember  2,  1988. 

3.  "Forced,  Unbounded  Shear  Flows",  with  F.  Browand  and  C.M.  Ho,  invited  talk,  CHAOS  ’87, 
Monterey,  CA,  January  15,  1987,  Nuclear  Physics:  B.  Proceedings  Supplement  2,  pp.  139- 
158,  1987. 

4.  The  Structure  of  the  Turbulent  Mixing  Layer",  F.  Browand,  Phvsica  18D.  pp.  135-148,  1986. 
(Invited  presentation.  International  Conference  on  Solitons  and  Coherent  Structures). 

5.  "An  Experimental  Study  of  Entrainment,  Mixing  and  Chemical  Reaction  in  a  Plane  Turbulent 
Shear  Layer",  with  T.  Maxworthy,  J.C.  Lasheras  and  J.S.  Cho,  Proceeding  of  the  20th 
International  Svmo.  on  Combustion.  August,  1984. 

6.  The  Origin  of  Three-Dimensional  Vortical  Structures  in  a  Plane,  Free-Shear  Layer",  with 
T,  Maxworthy,  J.C.  Lasheras  and  J.S.  Cho,  J.  Fluid  Mech..  1986. 

7.  The  Turbulent  Mixing  Layer  Geometry  of  Large  Vortices",  with  F.  Browand  and  T.R. 
Troutt,  J.  Fluid  Mech.  158.  pp.  489-509,  1985. 

8.  The  Mixing  Layer  An  Example  of  Quasi  Two-Dimensional  Turbulence",  with  F.  Browand 
and  C.M.  Ho,  Journal  de  Mecanioue  Theorioue  et  .Appliouee.  pp.  99-120,  Special  Issue  1983. 


Index  of  Reports 

VI.  Development  and  Implementation  of  Image  Processing  Algorithms 

1.  "Performance  Analysis  of  Automated  Image  Processing  and  Grid  Interpolation  Techniques 
for  Fluid  Flows",  with  E.J.M.  Rignot  and  G.R.  Spedding,  USC.AE  143.  1988. 

2.  "Application  of  Neural  Networks  to  Three-Dimensional  Particle  Tracking  in  Fluid  Flows". 
E.J.M.  Rignot,  USCAE  146.  1988. 


CONTROL  OF  MIXING  IN  FREE  SHEAR  LAYERS 
(Task  2) 

Principal  Investigator  Ho,  Chih-Ming 
ABSTRACT 


The  main  thrust  of  the  research  is  to  understand  the  physics  of  the  coherent  vortices  and  to  develop 
techniques  of  controlling  the  evolution  of  the  deterministic  structures.  Fundamental  understanding 
of  the  entrainment  and  the  small  scale  mixing  processes  in  free  shear  layers  have  been  obtained 
during  the  research  period  12/01/81  -  9/30/88. 

ACCOMPLISHMENTS 


I:  Control  of  2-D  Coherent  Structures 

The  coherent  structures  are  developed  from  the  Kelvin-Helmholtz  instability  and  therefore  are 
sensitive  to  the  initial  perturbations.  We  have  found  that  the  almagamation  pattern  of  the  vortices 
can  be  altered  dramatically  by  perturbing  the  shear  layer  at  the  subharmonics  of  the  most  amplified 
frequency  (Ref.  I- 1, 2,3,5,  and  6]. 

The  momentum  [Ref.  1-4]  and  the  mass  [Ref.  1-3,  S]  transfers  across  the  two  streams  of  a  free  shear 
layer  are  caused  mainly  by  the  unsteady  evolution  of  the  vortices.  Hence,  the  manipulation  of  the 
vortex  merging  patterns  has  important  implications  for  many  technical  applications. 

II:  Instrumentation 

The  split  film  probe  was  a  common  instrument  for  simultaneously  measuring  the  streamwise  and  the 
transverse  velocity  components.  When  it  was  used  in  our  laboratory,  we  could  never  obtain 
satisfactory  results.  After  a  careful  examination,  we  have  conclusively  shown  that  the  cross-talk 
between  the  two  films  ruined  the  phase  information  between  the  two  velocity  components. 
Therefore,  the  split  film  probe  is  not  a  usable  instrument  [Ref.  II- 1]. 

Ill:  Lock-on  of  Wake  Flow 

When  a  wake  was  forced  by  a  monochromic  disturbance,  the  wake  instability  waves  locked  to  the 
forcing  at  a  certain  frequency  band.  We  have  identified  the  ratio  between  the  forcing  frequency 
and  the  most  amplified  frequency  of  the  wake  as  the  important  parameter  for  characterizing  the  lock- 
on  band  [Ref.  III-l]. 

IV:  Subharmonlc  Resonance 

When  the  free  shear  layer  was  forced  by  a  single  frequency,  many  flow  quantities  could  be  predicated 
by  the  linear  instability  analyses.  When  the  flow  was  perturbed  by  the  fundamental  and  the 
subharmonic,  the  amplification  rate  of  the  subharmonic  was  found  to  be  a  function  of  the  phase  angle 
between  the  fundamental  and  the  subharmonic[Ref.  IV- 1].  This  phenomenon  supports  the  non¬ 
linear  subharmonic  resonance  mechanism  proposed  by  Kelly  [1967]. 

V:  Streamwise  Vortices 

The  vorticity  of  streamwise  vortices  was  found  to  be  in  the  same  order  of  the  maximum  mean  shear 
of  the  mixing  layer.  Furthermore,  the  wavelength  of  the  streamwise  vortices  doubled  after  the 
merging  of  the  spanwise  structures  took  place  [Ref.  V-1,2].  The  increase  of  the  length  scale  has 
been  confirmed  by  the  Navier-Stokes  numerical  simulation  [Ref.  V-3]. 


VI:  Phase  Dccorrelatioo 


Near  the  trailing  edge  of  the  splitter  plate,  the  phase  of  the  passing  coherent  structures  has  a  definite 
relationship  with  the  forcing  signal.  The  phase  jitter  was  significantly  increased  after  the  vortex 
merging  [Ref.  VI- 1,2].  The  phase  dynamics  and  the  small  scale  transition  were  considered  as 
possible  candidates  for  the  phase  decorrelation.  After  a  careful  examination,  we  have  identified 
the  background  noise  as  having  an  influence  on  the  growth  of  the  subharmonic.  Therefore,  we  can 
conclude  that  the  phase  decorrelation  starts  with  the  vortex  merging. 

VII;  Mixing  Transition 

In  an  originally  laminar  plane  mixing  layer,  we  examined  the  process  of  generating  the  random  small 
eddies  in  the  flow.  The  fine  eddies  were  first  detected  in  the  core  of  the  streamwise  vortices  at  the 
point  where  the  merging  of  the  spanwise  structures  occurred  [Ref.  VII- 1,  2,  3,  4,  S  and  6).  In  other 
words,  the  vortex  merging  not  only  transferred  the  energy  to  the  subharmonic  range,  but  also  shifted 
the  energy  to  the  high  frequency  end. 

VIII:  Sound  Generation 

A  numerical  code  of  the  compressible  plane  mixing  layer  was  used  to  study  the  production  of  sound 
waves  from  the  flow.  We  concentrated  specially  on  the  Mach  number  range  that  the  hydrodynamic 
wavelength  is  about  the  same  as  the  acoustic  wavelength.  In  these  cases,  the  pressure  field  can  be 
calculated  from  the  shear  layer  to  the  acoustic  far  field  [Ref.  VIII- 1 }.  The  results  show  the  radiation 
patterns,  the  Mach  number  dependence,  the  mechanism  of  sound  production,  etc.  Many  of  the 
unknown  aeroacoustic  phenomena  were  revealed  by  this  powerful  numerical  method. 
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Summary  of  Work  Accomplished 

I.  Inflectional  Velocity  Profiles  and  a  Turbulent  Production  Mechanism 

The  major  discovery  found  during  this  contract  was  that  turbulence  production 

within  transitional  and  turbulent  boundary  layers  is  associated  with  instantaneous 

inflectional  velocity  profiles.  This  is  a  different  and  more  distinct  mechanism  than 

has  been  proposed  before  and  suggests  some  new  methods  of  controlling  the  turbulence 
production  as  discussed  by  Blackwelder(i988)  and  Swearingen  and  Blackwelder(1989). 

Essentially  this  new  idea  for  turbulence  production  evolves  around  the 

instability  of  a  shear  layer  with  an  inflectional  velocity  profile.  It  has  been  known 
since  the  time  of  Rayleigh  that  these  profiles  are  unstable  to  small  disturbances  and 
have  an  extremely  rapid  growth  rate.  The  growth  of  a  disturbance  is  so  large  that  laminar 
wakes  and  Jets  are  almost  impossible  to  obtain  experimentally.  Instantaneous  velocity 
profiles  obtained  in  transitional  and  turbulent  boundary  layer  during  the  course  of  this 
research  revealed  that  these  flows  have  inflection  points  also,  although  they  do  not 

occur  continuously  in  time  nor  everywhere  in  space.  Examples  of  these  profiles  and 
details  of  the  production  mechanism  are  given  in  Blackwelder  and  Swearingen(1989) 

attached  as  appendix  A.  They  show  that  these  profiles  exist  for  sufficiently  long 
periods  of  time  that  they  can  be  considered  as  quasi-stationary  and  hence  the  results 

from  steady  state  stability  theory  are  applicable.  They  also  show  that  the  profiles 

occur  at  distances  sufficiently  removed  from  the  wall  that  the  viscous  effects  can  be 
ignored. 

One  of  the  more  surprising  aspects  of  this  mechanism  is  that  inflectional  profiles 

occur  in  both  the  normal  and  spanwise  directions  in  transitional  and  turbulent  layers. 
This  is  documented  in  figures  7  and  8  of  Appendix  A  for  a  transitional  layer  on  a  curved 
surface  and  in  figures  10  and  1 1  for  the  turbulent  case.  Heretofore,  the  spanwise 

direction  has  been  neglected  in  the  mechanisms  for  turbulent  production;  however  these 
new  results  suggest  that  the  local  inflectional  regions  should  be  important  and  must  be 

incorporated  into  any  successful  theory  and  model.  The  spanwise  structure  is  especially 
important  in  attempting  to  explain  the  success  of  riblets  in  drag  reduction;  i.e.  the 

riblets  may  act  to  slow  the  streamwise  vortices  and  thus  reduce  the  inflectional  aspect 

of  the  velocity  profile  which  in  turn  would  reduce  the  production  of  turbulence. 

II.  Large  Amplitude  Laminar  Disturbances 

Another  phase  of  this  research  during  the  last  seven  years  has  been  to  study  the 
breakdown  of  a  growing  laminar  disturbance  into  a  turbulent  patch  of  fluid.  The  earlier 
work  of  Blackwelder  et  al.(1984)  and  Tso  et  al.(1984)  concentrated  upon  examining  the 
initial  conditions  of  a  small  perturbation  within  a  laminar  boundary  layer.  The 
importance  of  the  spatial  structure  of  the  initial  disturbance  was  examined  by  using 
various  hole  diameters  and  non-circular  shaped  holes  for  the  injection  of  the  minute 

disturbance.  It  was  found  when  the  initial  disturbance  is  a  pulse,  that  growing  wave 
packets  were  produced  only  for  the  smallest  diameter  circular  holes;  e.g.  0.1  mm  in 
diameter.  Larger  holes  and  non-circular  holes  tended  to  produce  a  different  disturbance 
that  broke  down  into  turbulence  much  earlier.  In  all  cases  the  initial  break  down  into 
turbulence  occurred  at  about  y  -  0.25  and  appeared  upstream  of  the  breakdown  at  y  «  1.15 
observed  by  Caster  and  Grant(Proc.  Roy.  Soc.  Lond.  A 347.  1975).  After  the  turbulent 
regions  had  a  chance  to  develop  downstream,  a  classical  turbulent  spot  was  observed  in 


turbulence  occurred  at  about  y  >  Q.2S  and  appeared  upstream  of  the  breakdown  at  y  >  I.IS 
observed  by  Caster  and  Grant(Proc.  Roy.  &c.  Lond.  A347.  1975).  After  the  turbulent 
regions  had  a  chance  to  develop  downstream,  a  classical  turbulent  spot  was  observed  in 

all  cases. 

The  temporal  aspects  of  the  initial  disturbances  were  studied  by  Tso  and 

Blackweider(1985)  by  comparing  the  pulse  disturbance  described  above  with  a  tone  burst 
consisting  of  4  cycles  of  a  sinusoidal  wave.  At  low  amplitudes,  the  tone  burst 
disturbance  evolved  into  a  Caster  type  wave  packet  downstream.  With  a  large  amplitude, 
a  classical  turbulent  spot  was  generated  almost  immediately.  At  intermediate 

amplitudes,  a  wave  packet  was  observed  initially  downstream.  However  another 

perturbation  was  also  observed  which  had  a  larger  lengthy  scale  and  faster  propagation 

velocity.  Further  downstream,  this  larger  disturbance  grew  to  greater  amplitude  and 

moved  ahead  of  the  wave  packet  which  was  travelling  at  the  speed  of  the  Tollmien- 

Schlichting  waves.  This  larger  disturbance  broke  down  into  a  classical  turbulent  spot 

and  the  wave  packet  remained  behind  it  as  seen  by  Wygnanski  et  al.(JFM,^,  505,  1979). 
This  observation  explains  the  existence  of  the  trailing  wave  packets  and  illustrates  why 

the  packets  are  correlated  with  the  generation  mechanism  as  observed  by  Wygnanski  et  al. 
Tso  et  al.(1988)  have  continued  to  investigate  the  breakdown  of  these  different  type  of 

disturbances  into  turbulent  spots. 


III.  Overhang  Region  of  a  Turbulent  Spot 

It  had  been  conjectured  in  the  literature  that  the  overhang  (of  nose)  region  of 

a  turbulent  spot  is  very  important  in  the  growth  of  the  spot  because  it  contains 

turbulent  fluctuations  that  could  disturb  the  unstable  laminar  boundary  lying  beneath 
it.  These  disturbances  would  then  be  amplified  by  an  instability  mechanism  in  the 
laminar  region,  grow  and  breakdown  into  turbulence.  The  shape  of  the  overhang  region  and 
other  aspects  of  the  leading  edge  of  the  spot  were  investigated  by  Cutmark  and 

Blackwelder(1987).  The  distance  that  the  overhang  extended  upstream  from  the  turbulent 
interface  at  the  wall  increased  continuously  as  the  spot  grew  downstream.  It  appeared 
that  the  extent  of  the  overhang  was  a  function  of  the  previous  history  of  the  spot  and  not 
necessarily  a  function  of  the  Reynolds  number. 

The  mixing  characteristics  within  the  turbulent  spot  were  examined  by  Swearingen 
and  Blackwelder(1985)  by  generating  a  fully  developed  spot  and  passing  it  over  a  heat^ 
wall  region.  The  old  turbulence  was  colder  fluid  and  the  newer  turbulence  generated  in 

the  vicinity  of  the  wall  was  warm.  The  results  indicated  that  the  new  turbulence  is  very 
dynamic  and  immediately  dominates  the  flow  field.  However  a  sharp  demarcation  between 
the  heated  and  non-heated  regions  could  not  be  obtained  due  to  the  conductivity  of  the 
plate;  thus  this  line  of  research  was  not  continued. 
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INTRODUCTION 

Much  of  the  research  on  turbulent  boundary  layer  during  the  past  twenty  years  has 
concentrated  on  the  coherent  eddies  that  control  the  flow  field.  The  la^e  scale 
structure  of  the  outer  flow  field  dominates  the  flow  features  there  and  contols  the 
entrainment  of  irrotatlonal  fluUl  into  the  turbulent  regions.  Near  the  wall  the  flow 
field  is  dominated  by  a  senes  df  events.  coQectlvely  called  the  bursting  phenomenon, 
that  are  responsible  for  the  turbulence  production  and  the  drag.  This  series  of  events 
consists  of  cotmter-rotatlng  streamwlae  vortices,  low  speed  streaks  of  fluid,  lift  up  of 
that  fluid  off  the  wall,  an  oscillation  of  the  lifted  low  speed  streak  and  the  rapid  growth 
of  that  oscillation  and  its  break  down  Into  turbulence.  This  sequence  has  been 
described  in  detail  by  Kline  et  aL  (1967).  WlDmarth  (1975).  Cantwell  (1981)  and  others. 
The  main  thrust  of  the  present  work  is  to  esamine  the  osdilatlon  stage  and  suggest  a 
mechanism  by  which  this  occurs. 

STABlLrrr  THBORT  for  INFUCnONALVELOCXXT  propurs 

One  of  the  eariiest  successes  In  stability  theory  was  due  to  Raleigh  (1880)  who  showed 
that  inflectional  velocity  profiles  in  an  invlscld  fluid  are  unstable  to  small 
disturbances.  Viscosity  has  only  a  damping  influence  in  these  flow  fields  as  discussed 
by  Lin  (1955).  MlchaUce  (1965)  anatysed  the  hyperbolic  tangent  profile  for  a  mean 
parallel  flow  field  given  by 

U(0  --^ttanhl^A)! 

where  Uq  is  the  maximum  velocity  difference  across  the  shear  layer  and  2A  is  its 
vortlclty  thickness  as  illustrated  in  figure  1.  He  assumed  a  spatially  growing  two 
dimensional  disturbance  in  the  form  of  a  travelling  wave  given  by 

u'(x4,U» 

v'(x4.t)  a  ^(ye**®“'^ 

a  is  the  complex  wave  nunflier.  a  *  Of  iO}.  and  o>  is  a  real  frequency  where  x  and  ^  are  the 
streamwise  and  cross  stream  coordinates,  u’  and  v*  are  the  respective  velocity 
perturbation.  His  results  yielded  a  broad  range  of  growing  harmonic  disturbances  with 
the  most  amplified  wave  number  being  OjA  >  0.40.  The  corresporuling  wavelength  is 
a  15.7A  which  is  long  compared  to  the  shear  layer  thlclmeas. 

From  Near-Wal )  Turbulence,  a  Symposium  held  May  16-20,  1988  in 
hubrovn Ik,  Yugoslavia  published  by  Hemisphere  Pub.  Corp. 


FV.  1.  Sketch  of  an  mviadd]^  unstable  velocity  pira61e  m  the  (  direction. 


More  surprising  however  is  the  growth  rate  of  these  disturbances.  Of.  Mlchalke’s 

analysis  shows  that  the  disturbances  grow  like  where  the  maalnuim  growth  rate 
was  found  to  be  » *0.23.  Thus  while  the  wavellke  disturbance  travels  one 
wavelength  downstream.  x-Xq  *  it  is  amplified  by  a  factor  of  36  and  its  energy 
increases  more  than  1000  fold!  For  compaziaon,  the  maximum  ampllflcatlcn  of  the 
most  unstable  ToUmien-Schllchting  wave  is  typlcalty  1.1  (Le.  they  ercperlence  a  10% 
Increase  in  amplitude)  during  a  journey  of  one  wave  length  downstr^m  because  its 
growth  rate  Is  an  order  of  magnitude  smaller.  Hence  (Usturbances  associated  with 
inviacid  ir^^Uonal  velocity  profiks  are  amplified  much  more  than  those  associated 
with  viscous  mstablUtlea.  In  some  transitional  flow  fields,  the  initial  viscous 
instability  amplifies  an  disturbance  In  a  manner  that  produces  an  inflectional 
profile  as  discussed  by  Blackwelder  (1983).  Ihls  opens  the  door  to  another  type  of 
instability,  usually  called  a  secondary  or  higher  order  instability.  If  the  secondary 
instability  Is  an  inflectional  type,  the  disturbances  will  grow  much  faster  than  the 
initial  disturbance  and  will  dominate  the  transition  process. 

One  of  the  more  interesting  characteristics  of  the  inviscld  instability  is  that  its 
ampUfleatlon  is  strongly  dependent  upon  the  shear  rate  of  the  inflectional  profile.  To 
Illustrate  this,  note  that  the  maximum  arxxplification  over  the  distance  X2*xj  Is  given 

by 

eO.23  (X2-Xi)/A  (1 

Ihus  the  smaOer  the  length  scale  A.  the  larger  the  growth  over  the  same  diatance.  Figure 
1  shows  that  A  is  inversely  proportional  to  the  maximum  shear  of  the  profile.  Le. 

^  }h 
j  m  2A 

Eliminating  A  between  these  two  expresalons  yields  a  growth  of 


Consequently  the  maximum  amplificaUon  experienced  by  a  disturbance  over  the 
distance  X2'X^  downstream  depends  strongly  on  the  mean  shear  at  the  point  of 
inflection. 

The  above  arguments  following  Mlchalke's  analysis  have  been  for  a  steady  two- 
dimensional  l^ectlonal  profile  in  a  parallel  flow.  To  apply  these  arguments  in  a 
transitional  or  turbulent  flow,  the  restrictions  imposed  by  these  assumptions  must  be 
addressed.  The  effects  of  non-parallel  flow  have  been  examined  by  Caster  et  aL  (1985). 
They  found  that  when  the  local  value  of  the  mean  flow  was  used  for  the  stability 
calculations,  the  theory  accurately  predicted  the  measured  fluctuations  even  though  the 
flow  was  nonlinear.  Greenspan  and  Benney(1963)  have  studied  the  effects  of  an 

unsteady  mean  flow  by  allowing  the  thickness,  d.  to  vary  sinusoidally  with  time.  The 
forcing  ^riod  was  much  larger  than  the  time  scale  of  the  instability.  They  found  that 
the  unsteadiness  Increased  the  growth  rates  so  that  the  energy  id  the  disturbance 
increased  by  10.000  over  one  wavelength! 

Another  limitation  of  the  above  analysis  is  that  it  pertains  only  to  unbounded  flow 
domains.  However  Huerre  (1983)  has  considered  the  same  inflectional  velocity  profile 
confined  by  two  side  walls  located  at  ^  a  ±d.  He  assumed  the  normal  velocity  on  the 
walls  was  zero  but  allowed  a  slip  condition  to  occur  to  slzDpllfy  the  analysis.  He  found 
that  if  d  >1.2d.  the  linear  stability  characteristics  were  slmflar  to  the  inviscld  theoy 
and  the  wan  had  little  effect  on  the  growth  of  the  disturbances. 

The  two-dtmenaionaltty  of  the  mean  flow  has  not  yet  been  addressed  and  is  stffl  an 
open  questum.  However.  Nishioka  et  aL  (1980)  applied  the  invladd  two-dimensional 
theory  to  an  Inflectional  U(y)  velocity  prtdUe  in  a  transitional  channel  flow.  The 
Instantaneous  velocity  profile  was  formed  by  the  ToUmien-SchUchting  waves  and  was 
strongly  three-dimensional.  ExpertmentaUy  they  observed  that  a  secondary 
instability  appeared  at  the  location  of  the  inilectlon  point.  To  *««««"»  this  in  more 
detail,  they  applied  the  two-dhnenalonal  Invlacld  calculation  discussed  above  and 
found  that  it  predicted  the  measured  wavelength  and  ftigewfiinrtion  The  measured 
growth  factor  was  approximately  60%  of  the  theoretical  value.  Similar  results  have 
been  observed  in  transitional  boundary  layers  on  flat  and  concave  plates.  The  reduced 
growth  rate  may  be  due  to  the  spreadtog  energy  into  the  third  dimension  a^ich  is  not 
accounted  for  In  the  two-dlmenaioiial  theory.  This  effect  indicates  that  the  disturbance 
would  be  amplified  by  a  factor  of  9  after  travelling  one  wavelength  downstream  rather 
than  by  36  mentioned  earlier.  The  net  result  is  that  even  though  the  growth  rate  may  be 
reduc^  by  the  three-dtmenslonaUty.  It  Is  stiU  very  large. 

tOBrnMRmnjasurraiiAocmcM/EWMj, 

The  Gbitler  tnatablUty  is  one  member  of  a  class  of  centrifugal  instabilities  including 
the  Tqdar,  Dean.  etc.  instabilities.  These  instabilities  arise  as  a  result  of  an  imbalance 
between  a  centrifugal  force  imposed  on  the  fluid  and  the  restoring  pressure  force.  The 

Gbrtler  instability  develops  within  a  boundaty  layer  developixu  on  a  concave  wall  The 
curvature  of  the  wall  imposes  a  centrifugal  force  withm  the  fluid  that  acts  to  displace 
fluid  particles  toward  the  wall  The  pressure  force  acts  to  restore  the  particles  to  their 
orlgtnal  position.  This  class  of  mstabllttes  is  mviacid  m  nature  and  the  viacoatty  only 
acts  to  inhibit  the  motion.  Le.  provides  a  damping  force.  Just  as  in  the  inflectional 
profile  Instabilities.  Therefore  if  the  centrifugal  force  is  smaller  than  the  pressure 
force,  the  flow  is  stable.  However  if  the  cemnfugal  force  sufficiently  exceeds  the 
pressure  force  such  that  the  viscosity  can  no  kmger  damp  the  motion,  the  flow  field  is 
unstable  and  small  disturbances  win  be  amplified. 


Gbrtlet<1939)  Orst  analysed  this  problem  using  temporal  stability  theory  and  found 
that  the  governing  parameter  for  this  flow  was 


This  dimensionless  number  is  known  today  as  the  Gbrtler  number.  If  it  exceeds  a 
critical  value  of  the  order  of  unity,  the  flow  Is  unstable  and  develops  counter-rotating 
stieamwise  vortices  lying  withm  tte  boundary  layer.  These  vortices  are  periodic  in  the 

spanwtse  direction  and  are  characterized  by  a  qranwise  wave  length  Hall  (1983)  has 
shown  that  for  vortices  with  wavelengths  comparable  to  the  boundary  layer 
thickness,  the  growth  of  the  boundary  layer  cannot  be  neglected  in  the  analysis. 
Consequently  a  unique  neutral  stability  curve  does  not  exist  Imt  rather  depends  upon 
the  history  of  the  boundary  layer.  Nevertheless.  Hall  found  that  each  neutral  stability 

curve  had  a  minimum  at  Go  >  1.5. 

In  a  spatial  analysis  of  the  Gbrtler  Instability  m  a  parallel  flow.  Smith  (1955)  asaumed 
the  disturbances  on  a  concave  wall  grew  by  the  factor  as  they  travelled 

downstream  from  x^  to  X2.  3  is  the  growth  rate  and  varies  alowty  m  the  downstream 
distance.  As  the  fluid  moves  downstream  from  the  leading  edge,  the  disturbances  are 
initially  damped.  Le.  p  <  0.  They  crossed  the  neutral  stability  curve  at  P  «  0  and  ate 
unstable.  Le.  9  >  0.  further  downstream.  To  compare  this  growth  rate  with  the  lnviacld 
one.  9  ^  determined.  Liepmann  (1945)  has  shown  that  Gbrtler  instability 

transitions  to  turbulence  when  Go  ■  9.  The  value  of  ampllflcatton  rate  aUghtly 

before  transition.  Le.  Gb  •  7.  is  9^^  3  accordtng  to  Smith.  For  this  value  of  the 

Gbrtler  number.  Rcq  ranges  from  230  reported  by  Btypes  (1972)  to  360  from  Swearingen 
and  Blackwelder  (1987).  For  Reg  •  300  and  with  the  boundary  Ityer  thickness.  8. 

equal  to  7.56  as  in  a  Blaatua  layer,  the  arxqiliflcatlon  becomes  Thus  the 

streamwise  vortical  disturbance  wiD  double  its  amplitude  after  travelling  108 
downstream.  For  compariaon.  the  most  unstable  ToDmlen-Schlichting  wave  at  Reg  » 

300  has  a  wavelength  of  88  and  an  astyUficadon  fiKtor  of  These  waves 

are  amplified  approximately  80%  after  travelling  108  downstream.  Hence  at  this 

Reynolds  number,  the  Gbrtler  instability  has  a  slightly  greater  amplification  than 
Tniiinton  A»hiirhHng  waves.  The  coBdusion  IS  that  the  growth  rates  in  both  boundary 
Iqrer  flows  are  at  least  an  order  of  magnitude  less  than  that  of  the  inflectional  profile 
instability  diacussed  earlier 


s  ;  1  »,6^.  Kf  .'  JiOJ  4  '.f 


The  above  linear  theory  only  describes  the  initial  growth  of  a  small  disturbance  on  a 
concave  waU.  Like  most  linear  theories,  it  can  say  nothing  about  the  later  growth 
stages  and  breakdown  of  the  perturbaflon.  In  a  visual  study  m  the  vortices  on  a  curved 
plate  Bippes  (1972)  has  shown  that  the  perturbation  velocity  profiles  agreed  wefl  ^th 
the  computed  elgenftmctions  and  the  growth  rates  were  compaarable  to  the  values  given 
by  the  linear  theory.  Hfe  also  found  that  the  disturbances  ffiw  to  an  amplitude  of  a  few 
percent  of  the  streamwise  mean  velocity  vdule  maintaining  their  streaorwiae  vortical 
structure.  However,  after  a  limited  donain  of  linear  growth,  the  flow  field  devdoped  a 
strong  oadUatory  motion  before  breaking  down  mto  turbulence. 
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Fig.  2.  Instantaneous  visualizations  of  the  Goitler  instability  and  tranattkm  on 
a  concave  wall  The  concentrated  smoke  marks  the  km  speed  streaks 
and  the  smuous  motion  is  a  result  of  the  Inflectional  instability. 

This  aspect  of  the  motion  was  considered  to  be  important  in  the  transition  process 
because  it  represented  a  shift  fttsn  a  primary  to  a  secondary  instabiUty  m  the  flow. 
Since  ««wuiar  observations  have  been  made  in  the  wall  region  of  turbulent  flows  (see 
below),  an  experiment  was  set  up  to  measure  the  velocity  conponents  in  more  detail  to 
study  this  phenomenon.  The  speclaily  designed  wind  tunnel  used  for  this  experiment 
had  a  curved  test  section  with  a  cross  sectional  area  of  15  x  120  cm.  A  boundiiuy  layer 
was  developed  on  the  concave  wall  having  a  radius  of  ciuvature  of  R  «  3.2m  with  a  free 
stream  velocity  of  5.0  m/sec.  A  smoke  wire  was  used  to  visualize  the  flow  field.  The 
streamwlse  vdodty  was  measured  with  rakes  of  hot>wires  aligned  in  the  spanwlse  and 
normal  directions.  Thus  the  U(z.t)  and  U(y.t)  profiles  could  be  measured 
Instantaneous^.  The  ««gnai«  were  digitized  and  stored  on  disc  for  later  processing. 
Further  details  of  the  experimental  arrangement  may  be  foimd  in  Swemingen  and 
Blackwelder  (1987). 

At  a  freestream  velocity  of  5  m/sec.  the  critical  Gbrtler  number  occturred  at 
approadmately  10  cm  downstream  of  the  leading  edge  of  the  curved  plate.  When  the 
smoke  srtre  was  placed  into  the  lower  regions  of  the  boundary  layer  up^ream.  Le.  at  y  a 
0.258.  andx  a  20  cm,  the  resulting  flow  downstream  appeared  as  in  figure  2.  By  the  time 
the  was  convected  downstream,  it  had  also  been  coagulated  into  elongated 

regioru  as  seen  in  the  photo.  Swearingen  and  Blackwelder(1987)  showed  that  the 
regions  of  strong  smoke  concentrations  seen  In  the  figure  are  also  regions  of  low 
streamwlse  velocity.  They  also  showtKI  that  the  ensuing  dynamics  were  concentrated 
around  and  directly  relat^  to  the  low  speed  regions. 

A  sketch  of  the  streamwlse  vortices  and  the  low  speed  regions  is  shown  m  figure  3.  The 
instability  produces  the  counter-rotating  streamwlse  vortices  with  two  vortices 
per  spanwlse  wavelength.  The  strength  of  these  vortices  is  quite  small.  To  illustrate 
this,  consider  a  mean  parallel  flow  U(^  with  a  disturbance  flidd  u  such  that  the  total 
velomty  is 

u(x.y.z)  >  U(y)  ♦  u'(x,y.z) 

v(x.y.z) «  v'(x.y.z)/Re  (4 

w(x.y.d  •  w’(x.yjt)/Re 


where  Re  Is  the  RqrnoUi  number.  Ihleaceltiiglenresu'.T'.andw'aftbeenieoRlerae 
suggested  IqrOlpihiis  and  Stuait  (1973}  and  reified  hfnxyan  end  Sane  (1982).  When 
the  rerticity  oomponents  are  obtatned  ftem  these  eqiiattnns,  it  is  essdy  seen  that  the 
dlstuibanoe  streamwiae  rectlctty  Is  amsDer  than  both  ay  and  cBx  hy  the  factor  Re* 
Hence  eren  though  the  flow  field  is  generally  deacxihed  as  a  boundazy  hqrer 

with  embedded  growing  streamwise  reztices.  they  contain  ooly  small  amounts  of 
atreamwise  voztlclty. 

The  above  illustration  indicates  that  the  phase  of  coherent  eddies  la  much  more 
Important  than  the  amplitude  of  their  velocity  or  vmtlcity  fields.  It  also  su^ests  that 
detection  of  coherent  structures  should  not  be  baaed  on  amplitude  alone  without 

conaldeztag  phase  mfonnafion  also.  In  the  Gbitler  vortex  case,  the  vortices  do  generate 

dU  dU 

large  vertical  disturbances,  namely  oy  •  '^snd  «os  ■  But  aiiy  detection  scheme 
baaed  on  amplitude  alone  would  probably  £m  to  discern  the  underlying  eddy  structure, 
namely  the  streamwise  vortices.  Ofcourse  the  reason  the  vortices  are  so  efEectlre  is 

dU 

that  thqr  exist  In  a  region  of  strong  mean  shear.  and  thus  a  very  mull  ot^  rotatton 
can  produce  large  disturbances  acroas  the  spaiL 
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4.  Olustndon  of  the  relative  phase  of  the  dlatuxbance  velocity  components 
at  poattlons  belotr  the  centers  of  the  streannviae  vortices.  The  velocity 
amplitudes  depend  on  y  and  axe  not  to  scale. 

Between  the  vortices  in  Qgure  3  are  alternating  regions  of  low  axul  high  speed  Quid. 
These  regions  are  a  natural  result  of  the  cotmter>rotatlng  vortices  lying  la  a  strong 
Qow.  The  velocity  Qeld  associated  with  the  vortices  has  a  normal  velocity 
rwmpfttwn*  With  a  relative  maztxxnsn  or  minimum  between  the  vortices  at  poatflons  a 
and  c  respective^  as  seen  in  figure  4.  Simultaneously  the  fiuid  ths  centers 

of  the  vortices  has  a  spaxxwlae  velocity  component  that  alternates  m  sign  with  the 
vortices  so  that  its  modan  is  directed  towards  the  region  of  positive  normal  velocity. 
The  combined  action  of  two  vortices  is  to  diqiJace  lower  ^>eed  Quid  adgmaOy  near  the 
wall  Into  a  region  lying  between  the  vortices  shown  at  position  a  in  figure  4.  This 
creates  the  eloEgated  low  speed  regions  seen  In  figure  3.  On  the  other  hand,  the  higher 
speed  Quid  lying  above  the  centerline  of  the  vortices  experiences  a  relative  motion  that 
forms  concentrations  of  high  speed  Quid  between  alternate  vortices  at  poaltlona  c  m 
figure  4. 

One  of  the  principal  resulta  of  this  vortical  matum  and  the  uniipie  aspect  of  the  Gortler 
instability  la  that  it  creates  mfiectlonal  velocity  profiles  in  a  fiow  field  m  vHUch  there 
were  none  before  the  disturbance  started  to  grow.  It  is  interesting  to  note  that  a 
spanwise  inflectional  profile  wiU  develop  as  soon  as  the  vortices  begin  to  grow.  This  is 
best  tUustrated  by  considaing  the  usual  velocity  perturbatlona.  Le. 

u*  (x.y.d  >  *  ^  cos  oz  e^ 

V  (x.yjd  ■  v(y)cosazeP*  (5 

w*  (x.y.d  *  •  sin  os 

As  soon  as  this  disturbance  begins  to  grow,  there  are  inflection  points  in  the  u(z)  profile 
at  cs  a  ±(n  .  n  «  0.1.2...  The  strength  of  the  strain  rate  at  these  points  is 

I  ^1 

and  it  increases  in  magnitude  with  the  growth  of  the  instability.  Note  that  the  sign  of 
the  strain  rate  la  not  essential  in  the  inviacld  instability;  only  the  absolute  value  is 
Important.  Later  m  the  development  of  the  flow,  inflecuonal  profOes  develop  as 
prMllcted  and  obeetved  by  sevmal  tnvesdgatois.  These  inflectional  profiles  open  the 
door  to  an  inviadd  in^abiUty  as  discussed  above. 


Prom  top  to  bottom  are  the  Mo-contoun  of  a)  tbe  mean  vdoclty.  b)  Uie 
nns  Tdodty.  c)  the  nonnel  ahear  au/df  and  d)  ttie  apamviae  ahear  SU/aa. 
The  data  were  taken  at  7  a  00  cm  and  axe  aormaUxed  wtdi  the  ftee  stream 
velocttjr  and  the  eelocltf  gradient  at  the  wall 
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similar  tn  ft  jure  5  StX 


ng.  7.  Two  examples  of  the  spaxxwise  velocity  proOle  having  inllectlan  points 
and  laige  shear. 

local  dU/dz  maxima  occurs  at  the  identical  location  of  the  maxima  m  the  nns  eontours 
and  their  contour  shapes  are  very  •tmOar  (compare  with  figures  fib  and  fid).  Even  the 
new  local  maximum  near  the  wall  at  x*  100cm  in  the  fluct\iatlon  contours  is  also  found 
in  the  iso-shear  plot.  Hence  it  it  obvious  that  the  velocity  fluctuations  are  associated 
primarily  with  t^i.;  spanwise  shear  rate.  dU/dz. 

Examples  of  the  mean  spanwise  velocity  profiles  obtained  at  the  y  position  near  the 
maximum  in  the  dU/dz  contours  are  shown  In  figure  7  for  x  ■  80  and  100cm.  The  mean 
spanwise  averaged  velocity  is  approximately  O.TUw  in  both  cases  and  the  maximum 
shear  is  typically  0.5  dU/dyfw  Illustrate  the  stability  theory  predictions,  the 

theoretical  wavelength  associated  with  these  profiles  should  be  »  15. 7d.  From  figure 
7  at  100  cm  the  shear  Ityer  thickness  is  approximately  0.25cm.  This  yields  A,  3.9cm 
which  compares  favorabty  with  the  wavelength  seen  in  figure  2. 

Similar  profUes  in  the  normal  direction  are  shown  in  figure  8  at  four  downstream 
locations.  The  inflection  points  are  noted  and  it  is  easily  seen  that  the  inflections  move 
rapldty  away  from  the  wall  as  the  flow  moves  downstream.  The  velocity  at  the 
inflection  points  is  typically  0.65Uot. 

WALLSEGKJNS  OP  BOUNDED  TVRBOLENT  SHEAR  FLOWS 

The  waD  r^lon  of  bounded  turbulent  shear  flows  is  dominated  by  the  bursting  process 

which  is  composed  of  a  series  of  events  similar  to  that  described  above  for  the  Gbrtler 
problem.  Namely  there  are  streamwiae  vortices,  elongated  regions  of  low  speed  fluid 
called  low  speed  streaks,  lllt-up  and  oscillation  of  these  streaks  and  a  rapid 
disintegration  of  the  coherent  motion  followed  by  strong  mixing. 

The  role  of  the  streamwiae  vordces  in  the  turbulent  flow  is  less  dear  than  m  the  Gbrtler 
case.  Indeed,  there  is  no  genen^  accepted  definition  of  such  vortices  found  m  the 
literature.  Let  us  make  the  foUowtng  definition: 

A  vortex  consists  of  a  coherent  vortical  parcel  of  fluid  such  that 
the  instantaneous  streamlines  m  the  plrae  perpendicular  to  the 
vortex  lines  are  closed. 


By  the  time  the  flow  has  reached  100cm  downstream  in  the  experiment,  the  velocity 
field  was  strongly  nonlinear  as  evidenced  by  the  variations  in  the  measured  streamwlse 
velocity  component.  Figures  5  &  6  illustrate  the  iso-veloclty  contours  measured  over 
two  wavelengths  in  the  spanwtse  direction  at  two  downstream  locations,  x  s  80  and  100 
cm.  Two  low  speed  regions  at  z  2  -0.8  and  ■♦•0.8cm  are  evidenced  by  the  relative  lower 
speed  fluid  lyl^  at  the  elevated  regions.  The  U(y)  s  0.9  contour  lies  well  beyond  the 
nominal  boundary  layer  thickness  defined  by  U(5^)  »  0.99U«.  The  growth  of  the  low 
speed  region  has  been  quite  large:  the  U  »  0.9Uw  contour  has  been  displaced  outward 
from  1.25  to  1.75  where  is  the  corresponding  Blaslus  layer  thickness  at  each 
location.  Hence  the  disturbance  has  a  large  growth  rate  compared  to  the  nominal 
laminar  boimdary  layer.  The  displacement  thickness  at  these  locations  was  measured 
to  be  approximately  8*  *  0.7cm  compared  with  the  Blaslus  value  of  5*  =  0.28cm.  The 
corresponding  higher  speed  regions  at  z  z  0  and  ±  1.6cm  have  thixmed  the  boundary 
layer  considerably  such  that  8*  -  0.1cm  there.  In  these  locations  the  shear  is  especially 
strong  near  the  compared  with  the  shear  in  the  lower  speed  r^ions. 

As  noted  earlier,  the  Gbrtler  instability  grows  spatially  downstream  and  has  no 
temporal  dependence.  However  oscillations  do  develop  downstream  as  seen  tn  figure  2 
which  indicate  the  presence  of  a  secondary  instabfilty.  The  temporal  fiuctuations  were 
measured  and  their  iso-contours  are  given  in  figures  5b  and  ^  at  the  two  locations. 
Note  that  the  maximum  at  each  position  occurs  on  the  sides  of  the  low  speed  streak  and 
not  near  the  top.  The  maximum  magnitude  increases  from  roughly  0.03U«  to  0.08U« 
bewteen  80  and  100  cm  downstream.  The  total  energy  associated  with  the  temporal 
oscillations  is  proportional  to  the  square  of  the  u*  values  integrated  over  the  y-z  plane. 
Clea^  this  energy  has  increased  more  than  an  order  of  magnitude. 

It  is  also  significant  that  the  lodl  of  the  lao-contours  have  moved  outward  from  the 
waU.  The  maximum  at  x  a  80cm  is  at  O.Scm  and  the  corresponding  at  z  « 

100cm  is  at  1.15cm.  The  corresponding  Blaslus  thickness.  8^,  haa  only  increased  from 

0.78cm  to  0.86cm.  This  outward  movement  corresponds  to  an  average  normal  velocity 
greater  than  0.03Um  which  is  more  than  an  order  of  magnitude  greater  than  the  roa««ii« 
velocity  nonnal  to  the  plate. 

Another  important  feature  of  figure  6b  Is  that  as  the  primary  fluctuations  have  moved 
outward,  a  new  local  maximum  has  evolved  m  the  waU  region  at  z  ■  100cm.  At  llOan 
downstream,  this  maximum  was  apprazixxxate^r  0. 12UU  whereas  the  outer  maztmnm 
had  remained  constant  at  0.08U«».  After  110cm.  the  endre  layer  became  turbulent  and 
the  local  maxima  were  no  longer  dlacemable. 

The  reason  for  the  occuirence  of  two  maximum  nns  values  per  wave  length  is  apparent 
when  the  iso-ahear  contouza  are  exammrd  in  figures  5cd  and  6cd.  The  dU/dy  and  dU/dz 
gradient  fields  were  calculated  by  fitting  the  averaged  velocity  daU  with  a  spline  fit  and 
then  differentiating.  They  are  normalized  with  the  average  gradient  at  the  waU  which 
was  1680  and  2300  sec'^  at  x  «  80  and  100cm  respectively.  The  dU/dy  iso-shear 
oontouia  are  presented  in  figures  5c  and  6c  arxl  the  corresponding  dU/dz  contours  in  5d 
and  fid  at  the  two  streamwlse  locations.  The  nonnal  shear  dU/^  has  its  at 

the  waU  as  expected,  but  also  develops  a  strong  local  maxima  above  the  low  speed 
streak.  This  maxima  moves  oirtward  as  it  travels  downstream  and  is  reduc^  subtly 
in  magnitude.  Distinct  regions  of  negative  shear  are  seen  at  x  >  100cm  due  to  the  low 
speed  flidd  riding  above  hl^er  velocity  fluid  as  seen  m  figure  6a. 

The  spanwlae  shear  rate.  dU/Sz,  has  a  rfiatiTirt  minimum  atvi  maximum  on  both  sides  of 
the  low  speed  streaks.  Symmetry  indicates  the  minimum  and  maximum  must  have  the 
same  majpntudes:  howew  for  these  two  streaks  studied,  the  negative  amplitude  was 
alwiqrs  slightly  less.  The  maximum  mamiitude  of  3U/dz  is  comparable  to  that  of  dU/dy 
at  X  a  80cm  and  exceeds  dU/dy  at  x  «  IOOczxl  The  loc&i  maxima  of  the  spanwlae  shear 

2 

Increases  from  0.4  to  0.5u^/v  downstream  In  contrast  to  the  decrease  in  dU/dy.  The 
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a  Evolution  of  the  dU/^  inflectional  velocity  profiles  in  the  x  direction 
along  the  centerline  <I  a  low  speed  region. 

This  definition  emphasizes  the  coherent  aspects  of  ihe  motion  which  are  associated 
with  and  due  to  Its  phase  and  not  necessarily  due  to  its  amplitude.  Thus  it  includes  the 

weak  Gbrtler  streamwise  vortices,  vortex  rings,  trailing  wing  tip  vortices,  etc.,  but 
excludes  the  mean  vortex  lines  in  a  two-dimensional  shear  flow. 

Evidence  of  such  vortices  in  the  wall  region  has  been  assumed  by  many  authors, 
however  very  little  quandatlve  ita  exists  supporting  their  exlstance.  Most  at  the  data 
that  are  available  are  visual  results  such  as  those  m  Kline  et  aL  (1967)  and  Smith  and 
Schwartz  (1983).  Both  of  these  authors  have  used  hydrogen  bubbles  to  visualize  the  flow 
field  and  have  shown  that  often  the  hydrogen  bubbles  are  coagulated  into  a  confined 
region  with  definite  rotation  perpendlciilar  to  its  principal  a:ns  consistent  with  the 
alK>ve  definition  of  a  vortex.  Smith  and  Schwartz  specidated  that  the  observed 
streamwise  rotations  were  due  to  the  legs  of  horseshoe  type  vortices.  Acalar  and  Smith 
(1987ab)  have  set  up  two  model  laminar  flow  fields  to  study  the  vortical  nature  of  the 
motion  near  a  aoUd  boundary.  In  the  first,  they  studied  the  flow  over  a  hemisphere  and 
the  second  was  the  flow  over  lower  speed  fluid  injected  throu^  a  streauwise  slit  in  the 
wan.  In  both  cases  they  observed  the  periodic  formation  of  repeatable  horseshoe 
vortices.  The  structure  was  studied  by  placing  the  hydrogen  bubbles  at  different 
positions  and  photographing  the  vortices  at  different  times  after  their  formation. 
Their  visual  results  are  very  ***"*!*«•  to  photographs  obtained  in  the  turbulent  wafl 
region  and  surest  that  similar  vortices  appear  there  also.  Kim  and  Moin  (1986)  have 
observed  similar  type  structiues  in  their  numerical  data  by  tracing  vortex  lines.  By 
examining  regions  of  strong  streamwise  vortlcity.  they  conduded  that  the  associated 
vortex  lines  rarely  remained  paraUel  to  the  waU  for  distances  greater  than  100  v/u^. 
This  restilt  is  consistent  with  the  above  observations  by  Smith  and  Schwartz. 

The  conditionally  averaged  data  of  Blackwelder  and  Eckehnarm  (1979)  and  of  Kim 
(1983)  also  support  the  Idea  that  streamwise  vortices  exist  in  the  waU  region.  Even 
thou^  their  results  yielded  symmetrical  vortices,  the  symmetry  was  a  result  of  the 
techmque  and  does  not  imply  that  symmetrical  pairs  win  occur.  Guezennec  et  aL  (1987) 
have  shown  that  the  vortices  rarely  appear  in  a  symmetrical  maxmer.  but  one  vortex  is 
usually  much  larger  than  the  other.  Obviously  spanwise  symmetry  requires  that  the 
average  streamwise  vortlcity  be  zero. 
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F^.  9.  Visualization  of  the  low  speed  streaks  and  eddy  structure  m  the  wall 
r^n  made  visible  by  a  hydrogen  bubble  wire  at  y*** «  12. 

The  most  ubiquitous  structures  of  the  waU  layers  are  the  low  speed  streak.  They  hsnre  a 
width  of  20v/ut.  an  average  spanwise  spacing  of  lOOv/u^.  and  have  steamwise  lengths  up 
to  and  exceeding  1000v/u.t.  The  low  speed  streaks  are  most  readily  observed  by 
hydrogen  bubbles  as  seen  in  figure  9.  The  hydrogen  bubble  wire  Is  seen  at  the  top  left 
hand  comer  of  the  photograph  and  was  pulsed  at  40  Hz:  The  whe  was  at  y-f  >  12  and  the 
mean  flow  velocity  was  10  cm/ sec.  The  low  speed  streaks  are  readily  recognizable  by 
the  closeness  of  the  bubble  lines  and  the  higher  concentration  of  bubbles.  The  spanwise 

scale  of  the  low  speed  regions  is  roughs  20v/ui.  The  spanwise  scale  of  the  higher  speed 
fluid  is  much  larger,  typically  about  SOv/u^,  so  that  the  spanwise  spacing  between  low 

speed  streaks  is  typically  SOv/u,.  All  of  these  dimensions  are  random  variables  and 
hence  only  average  values  can  be  used  for  the  scales. 

Since  the  pulsed  line  of  bubbles  is  a  time  line,  its  z  posiUon  downstream  of  the  bubble 
wire  is  proportional  to  the  Langranglan  streamwlse  velocity  component.  A  hi^- 
Ughted  time  line  is  shown  in  figi^  9  to  illustrate  the  U(z)  velocity  profile.  Neglecmog 
the  difference  between  the  Lai^rangian  and  Euleilan  velocities  for  the  moment,  one 
can  see  that  the  U(z)  profile  has  a  strong  variation  across  the  field:  its  approximate 
minimum  la  0.2Um  and  its  corresponding  maximum  is  0.8U...  Also  of  importance  is 
the  fimt  that  due  to  this  variation,  there  are  points  of  inflection  in  the  profile  on  both 
sides  of  the  low  speed  streaks.  These  regions  are  possible  sites  for  an  inviscld 
tnflecdonal  instability  as  discussed  earlier.  In  addition,  the  shear  rate  there  is  found  to 
be  extremely  large.  Le.  of  the  same  order  of  magnitude  as  the  shear  rate  at  the  wall 
2 

which  is  u^/v. 

Ihls  phenomenon  was  studied  in  more  detail  in  a  turbulent  boundary  layer  at 
Reg  ■  2200  by  using  multiple  hot-wire  probes  spanning  the  y  and  z  directions.  Since 
th^  are  sensitive  to  the  streamwlse  velocity,  they  measure  the  U(^  and  U(d  profiles  as 
seen  m  figure  10.  Ihe  data  from  each  sensor  is  plotted  and  the  solid  line  is  a  spline  fit  of 
the  data.  The  spline  fit  was  differentiated  to  determine  the  locations  of  the  inflection 
points  which  are  denoted.  A  comparison  of  figure  9  and  10  indicates  that  the  shear 
rates  are  slightly  higher  in  the  visualization  results  probably  due  to  the  finite 
resolution  of  the  hot-wire  sensors.  Nevertheless,  it  is  obvious  that  there  are  many 
inflection  points  in  both  U(y)  and  U(z}  in  the  wall  region. 


F«.  10.  Instantaneous  turbulent  boundary  layer  profiles  fiom  hot-wire  seoaots 
of  the  streamwise  velod^  component  m  the  (a)  nonnal  and  (b)  span-wise 
directions.  The  data  potnta  are  indicated  by  (4)  and  a  spline  fit  yielded 
the  lines.  The  points  of  tnllectlon  (  4  )  were  obtained  orom  the  spline 

2  2 

function.  The  time  between  profiles  is  1.5v/u  m  ta)  and  l.(V/u  in  (b). 


Fi^ie  1 1  gives  another  example  of  the  ubiquitous  nature  of  the  inflectional  profiles. 
These  lesutts  were  obtained  fiom  the  numerical  stnnilatlon  of  a  turbulent  chas^  flow 
at  Re9  •  200.  The  data  were  searched  for  points  of  inflection  which  are  denoted  on  the 
figure.  Comparing  figures  10  and  11  suggests  that  the  lower  Reynolds  number  channel 
flow  has  more  Irdlectlon  points;  however,  this  is  most  probably  an  illusion  due  to 
different  techniques  used  for  the  two  data  sets.  It  is  quite  evident  that  both  the 
experimental  and  numerical  data  indicate  that  points  of  inflection  are  very  evident  in 
both  the  and  Utd  data. 

A  comparlaon  d  figures  lib  and  11c  indicate  that  inflection  pomts  are  much  more 
predominate  m  the  wan  reglan  at  y***  •  15  than  in  the  logarithmic  region  at  y'*’ «  66.  This 
agrees  with  the  earlier  idea  that  the  low  speed  streaks  are  concomitant  with 
inflectional  profiles.  Since  the  low  speed  streaks  do  not  extend  outward  beyond  y***  -  25. 
one  expects  that  fewer  Inflection  pomts  wiU  be  observed  above  that  level 

In  addition  to  the  inflectional  character  of  the  velocity,  several  other  coiuUtlons  are 
necessary  before  the  simple  stability  anaylsis  can  be  applied.  First  the  flow  shotild  be 
steady;  le.  the  changes  m  time  m  the  base  flow  must  have  a  much  longer  time  scale  that 
the  time  scale  of  the  instability.  Smith  and  Metzler  (1983)  h^  measured  the 
persistence  of  the  low  speed  streaks  and  find  that  they  have  an  average  life  time  of 
2  2 
480v/u^  with  some  streaks  persisting  up  to  2580v/u^.  These  values  are  much  longer 

2 

than  the  tune  scale  of  the  inatabiltty  vdilch  is  of  the  order  v/u^ .  Thus,  a  quaai-steady 


ng.  11.  Instantaneous  velocity  profiles  in  the  wall  region  from  the  direct 
numerical  simiilatlon  of  a  turbulent  ^annel  flow  m  (a)  the  normal 
direction,  (bl  In  the  spanwise  direction  at  y'*’ «  15.  and  (c)  m  the  qianwise 
dliectlon  at  y*^  s  66.  Ihe  mfiectlon  polixts  are  denoted  by  4  . 

state  exists  for  the  inflectional  instability.  Secondly  the  site  of  the  instability  should 
be  fhr  removed  from  a  solid  boundary  to  insure  Its  inviscid  character.  Huerre  (1983) 
has  shown  that  the  characteristics  of  the  instability  are  hardly  altered  by  the  presence 

of  a  wan  as  long  as  the  inflectional  region  is  removed  at  least  1.2A  from  the  waU.  Since 

2 

the  length  scale  of  the  inflection  is  5-lOv/u^.  Huerre’s  crtteria  appears  to  be  satisfied  for 

the  inflections  above  y*^  >  10.  The  third  condition  necessary  to  apply  the  stability 
criteria  is  that  of  two-dhnensionality.  No  theory  is  presently  available  to  ejqplaln  how 
win  effect  the  instability,  but  Nlshloka  et  al's  (1980)  results  discussed  earlier 
suggest  that  the  two-dimensional  theory  is  valid  m  highly  three-dimensional  fields. 

To  determine  the  growth  rates  of  the  incipient  instability,  the  data  were  explored 
further  to  obtain  the  strength  of  the  shear  layers  associated  with  the  inflectional 
prt^le.  Ihla  was  accomplished  by  first  locating  the  position  of  the  inflection  points 
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Fig.  12.  Conditional  probabilit]r  distribution  of  the  Telodtj  shear  at  the 
toflection  points  &xni  experimental  data  with  Reg  ■  2500.  a)  velocity 
shear  pcniendicular  totfaewallat  10<y^<2a  b)  Velocity  shear  pnaDd 
to  the  wan  at  y^  •  IS. 

withm  a  given  spatial  domain.  AU  of  the  inflection  points  were  included  in  the 
averaging  process,  although  some  could  have  been  escluded  using  Fjgrtaft's 
The  value  of  the  shear  rate  was  then  measured  at  each  of  these  points.  Thus  the 
conditional  probability  of  the  shear  rate  cotild  be  obtained  by  MtutpHng  the  at  the 

points  of  milectlan.  The  results  are  shown  in  flgures  12  and  13  where  the  shear  rates  are 
nonnahaed  by  the  mean  wall  values.  The  first  of  these  figures  is  obtained  firom  the 
experimental  hot«wire  data  base  used  in  figure  10.  The  normal  shear  dlstnbutlon  in 
figure  12a  is  asymmetrical  because  of  the  mean  shear  in  that  direction.  The 
probable  value  of  the  shear  is  approximate  zero  and  its  mean  value  equals  dU/dy  at 
that  location.  It  is  also  interesting  to  note  that  the  negative  gradients  are  often 
observed  as  reported  earlier  by  WlUmarth  and  Shazma  (1984).  They  also  observed 
dU/dy  gradients  which  exceeded  the  mean  shear  at  the  walL  The 
spanwise  shear  in  figure  12b  is  symmetrical  as  required  and  has  dual  p»ain«  at  a 
nondimenaional  shear  of  approximately  ±0.1. 

Numerical  data  obtained  from  Spalart  (1987)  were  also  analyaed  to  compare  with  the 
experimental  results.  The  same  processing  was  used  to  obtain  figure  13.  The  I^nolds 
number  of  this  data  is  Reg  •  670.  The  results  are  very  to  that  in  figure  12;  the 

mam  difference  being  that  the  magnitudes  of  the  probabilities  are  slightly  lower  and 
the  distribution  is  allghtty  wider  than  m  the  experimental  case.  This  dltlference  may  be 
due  to  the  spatial  averaging  imposed  by  the  hot-wires  and  the  coarser  grid  spacing  m  the 
experimental  case. 


n*.  13.  Sfcnilar  muttakW  in  FtgUK  13  but  ibr  the  d&ect  numerical  stmuletloo  of 
a  turbulent  hmindwy  li^at  Hag  *670. 


If  the  ntoAlea  are  indeed  unatable.  then  one  ahould  be  able  to  predict  the 

wavdengtbQftheanqdtaedoocinatlan.  Mlchalke’a  resulta  surest  that  the  wavelength 
win  be  15.7A.  Close  examination  of  figures  10  and  11  and  other  similar  data  indicate 
A  la  roughly  10  ±6v/u^  Thus  the  mfiectlonal  profiles  should  produce  disturbances 
with  wavelengths  between  75  and  225  v/u,.  Some  distinct  data  exist  in  the  literature 
which  support  t*»i*  prediction.  Dlnkelatdier  et  aL  (1977)  measured  the  pressure 
fluctuations  on  the  wan  and  found  wavellhe  patterns  propogatlng  downstream  with 
wavelengths  of 200v/Uf  traveUl^  at  typlcafly  12  -  15u^.  Wine  et  aL  (1967)  also  reported 

observations  of  wavelengths  of  240v/u,.  One  reason  there  are  few  data  of  this  type 
reported  in  the  literature  is  that  the  oaclUatlons  are  occurring  within  a  random 
background  and  thus  they  are  dlfilcult  to  distinguish  from  the  uncorrelated 
In  since  the  background  provides  the  initial  disturbance  for  the 

instability  and  the  energy  of  the  disturbance  can  grow  by  a  factor  of  ICXW  while 
travelling  one  wavelength  downstream,  an  identifiable  perturbation  can  probably  only 
be  observed  for  one  or  at  most  two  wavelengths  before  it  is  distorted  beyond  the  point  of 
recognition. 


DlSCUaSIOIf  AND  CONCLITSIONS 


The  pnmaxy  result  of  this  research  is  that  inflectional  velocity  profiles  are  a  prevalent 
and  ubiquitous  feature  of  the  wall  region  of  turbulent  flow  fields.  The  presence  of  low 
speed  streaks  in  this  region  has  been  reported  by  numerous  authors  and  is  recognized  as 
one  of  the  clearest  indlrations  of  the  existence  of  wall  layer  eddy  structures.  The  mere 
presence  of  low  speed  streaks  is  sufficient  to  guarantee  that  there  will  be  inflectional 
velocity  profiles  in  the  spanwlse  direction  because  U(z)  must  vary  from  high  to  low 
values  at^t  a  zero  mean  as  seen  in  figures  9  through  11. 

The  instantaneous  data  also  shows  that  there  are  numerous  inflection  points  in  the 
direction  normal  to  the  wall  as  well.  This  is  somewhat  surprising  because  of  the 
restraining  influence  of  the  wall,  but  has  been  verified  by  diflerent  authors  in  several 

flows.  The  density  of  the  d^U/dy^  s  0  points  is  comparable  to  the  d^U/dz^  s  0  points  in 
the  region  10  <  y"*"  <  20. 

It  is  interesting  to  note  that  an  '’inflection  point”  is  actually  an  intersection  of  a  locus  of 
such  points.  In  the  simplified  two-dimensional  tanh(y/A)  profile,  the  locus  of  the 
inflection  points  is  a  plane  at  y  a  O:  i.e.  an  inflectional  plane.  'The  resulting  instability 
waves  have  crests  parallel  to  the  inflectional  plane  with  their  crests  aligned  along  the 
third  axis.  l.e.  perpendlciilar  to  both  the  mean  flow  and  to  the  direction  of  the  gradient.. 
These  waves  are  convected  along  the  x  axis.  For  the  spatial  growth  case,  the  wave 
amplitude  increases  in  x  at  each  point  in  time. 

In  the  three-dimensional  cases  considered  here,  the  locus  of  inflection  points  will  be  a 
non-planar  surface.  For  oeample.  a  wake  has  a  locus  of  inflection  points  which  farms  a 
cylindrical  surface  aligned  in  the  mean  flow  direction  with  equal  amplitude  of  the 
shear.  dU/d^.  at  every  streamwlse  position.  In  fully  developed  turbulent  flow,  the 
inflectional  surface  is  much  more  convoluted.  For  example,  the  side  of  a  low  speed 
streak  win  always  be  such  a  surface  because  d^U/dz^  s  o  there.  If  the  streak  has  an 

Inflectional  profile,  the  inflectional  surface  mqy  extend  from  one  side  of  the  low 
speed  streak  over  its  top  to  the  other  side  engulfing  the  entire  low  speed  regtoa  If  the 

strength  of  the  shear  dU/d^.  la  constant  over  the  entire  inflectional  sheet,  each 
instability  wave  wifi  have  a  similar  growth  rate.  'The  waves  would  be  locally  parafiel  to 
the  inflectional  sheet  and  hence  would  have  the  shape  of  a  horseshoe.  'The  vortices 
generated  by  Acariar  and  Smith  (1987ab)  clearly  illustrate  this  phenomenon  in  a  well 
controlled  laminar  flow.  Similar  vorttees  were  often  observed  in  the  transitional  flow 
of  Swearingen  and  Blackwelder  (1987). 

The  Inflectional  surface  is  not  necessarily  of  infinite  extent.  For  example,  a 
momentary  disturbance  at  a  point  in  a  laminar  boundary  layer  must  have  an 
inflectional  surface  of  finite  extent  because  the  flow  at  la^e  distances  from  the 
disturbance  may  be  devoid  of  infiActlon  points.  In  tracing  out  the  inflectional  surfaces 
In  data  such  as  m  figures  10  and  11.  one  often  sees  that  the  suifime  ends  abruptly  at 
some  point  in  the  flow  field  althm^  it  may  continue  in  the  third  dimension.  In  the 
turbulent  case,  the  mfWrMnnai  surfaces  are  random  m  space  and  time  as  well 

The  locll  of  Inflection  points  are  important  because  they  present  a  mechanism  by 
which  turbulence  is  produced.  Naniely  the  shear  associated  with  the  inflections  is 
generally  quite  strong  which  suggests  that  an  Invlscid  instability  dominates  the  flow 
field.  ‘Through  this  mechanism,  the  unstable  inflectional  profiles  can  extract  energy 
from  the  instantaneous  base  flow  and  produce  a  rapidly  growing  oscillation.  ‘The 
wavelength  of  the  oscillation  scales  with  the  width  of  the  inflectional  instability. 
Since  the  scale  of  the  is  random,  so  is  the  scale  of  the  growing  disturbance: 

however,  the  scale  of  all  of  the  observed  oscillations  was  well  within  the  spectral 
domain  of  the  turbulence.  The  instability  grows  extremely  rapidly.  i.e.  the  two- 
dimensional  theory  implies  that  the  oscHlations  Increase  their  amplitude  by  a  factor  of 
30  (and  Its  energy  by  1(X))  as  it  travels  one  wavelength  downstream.  The  actual  growth 
Is  probably  tempered  by  energy  spreading  in  the  third  dimension  as  weQ. 


A  model  flow  fleld  emulating  the  turbulent  wall  region  was  studied  which  consisted  of 
Gbrtler  streamwise  vortices  in  a  laminar  boundary.  This  flow  had  the  advantage  that  it 
was  steady  and  thus  provided  a  better  teat  bed  for  studying  the  velocity  profiles  and 
oscillations.  Similar  low  speed  streaks  and  inflectional  profiles  were  observed  in  the 
model  flow.  Inflectional  profiles  in  the  normal  and  spanwlse  directions  yielded  regular 
and  repeatable  oscillations  with  wavelengths  consistent  with  the  predictions  of 
Mlchalke  (1965).  The  oscillations  grew  rapidly  in  the  streamwise  direction  and  broke 
down  into  turbulence  after  two  to  four  wavelengths.  In  this  flow  fleld.  it  was  clear  that 
virtually  all  of  the  turbulence  was  produced  by  the  mviscid  instabill^.  Not  only  were 
the  dyriamlcs  similar  between  the  model  and  tiirbulent  flows  but  the  length  and  velocity 
scales  were  also  comparable  when  normalized  with  the  friction  velocity  and  the 
kinematic  viscosity,  as  noted  by  Blackwelder<1983).  This  scaling  is  common  in 
analyzing  turbulent  wall  flows  but  has  not  been  used  frequently  in  transitional  flows. 
The  present  results  suggest  that  this  scaling  provides  a  means  of  comparing  these  two 
flow  regimes  in  more  detail  than  heretofore  and  indicates  that  their  dynamics  are 
Indeed  similar. 

Finally,  the  evidence  suggests  that  the  instability  may  be  more  successftjl  in  producing 
u  and  w  fluctuations  than  in  producing  u  and  v  oscillations.  Note  that  the  instability 
will  always  produce  an  oscillation  in  the  direction  of  the  base  flow.  u.  In  a  two- 
dimensional  case,  the  second  velocity  component  produced  by  the  instability  is  the 
component  parallel  to  the  direction  of  the  gradient.  Since  an  iriflection  point  is  alwai^s 
present  on  the  sides  of  the  low  speed  streaks,  one  eaqiects  these  to  produce  oscillations  of 
u  and  w.  Oscillations  involving  v  may  be  less  preralent  because  inflections  do  not 
always  occur  above  the  low  speed  streak.  However,  it  appears  t^t  by  the  time  that  large 
growth  is  observed,  the  locus  of  inflection  points  exteiids  upward  from  the  sides  a^ 
over  the  low  speed  streak.  In  this  case  the  observed  oscillations  wm  depend  upon  the 
integrated  growth  rates  that  the  disturbances  have  experienced  which  may  fiirvor  the 
spanwlse  raciUatlon  since  they  have  existed  over  a  longer  time  period.  In  addition,  the 
^wth  of  the  V  oscillation  w^  be  hindered  at  some  point  1^  the  presence  of  the  wall, 
thus  suggesting  that  the  spanwlse  oscillations  will  dominate  and  produce  a  larger 
oactllation  amplitude.  This  la  of  course  consistent  with  the  observations  that  w'  >  v*  in 
the  wall  xegum  of  turbulent  shear  flows. 

AU  of  the  above  are  devoid  of  Reynolds  averaging  and  its  implications.  In  particular, 
these  arguments  view  turbulent  production  as  an  instantaneous  event  n^ted  to  a 
particular  mechanism,  e.g.  the  xmstable  inflectional  velocity  profile.  The  reported 
observations  show  that  the  growing  oscillations  are  three  dimensional  and  produce  u. 
v,  and  w  fluctuations  in  a  random  manner  depending  upon  the  orientation  of  the 
inflectional  profile.  This  observation  is  not  inconsistent  with  Reynolds  averaging  but 
does  provide  a  dlflerent  description  at  the  production  process.  For  example  in  the  above 
formulation,  the  turbulent  energy  is  fed  directly  into  the  u.  v.  and  w  fluctuations.  The 
Reynolds  averaged  equations  show  that  the  average  turbulent  energy  is  only  fed  into  the 
u  fluctuations  and  then  la  redistributed  by  the  pressiue  into  the  other  components.  The 
difference  is  that  the  averaged  equations  do  not  allow  for  instantaneous  variations  in 
the  base  flow.  Le.  the  mean  averaged  flow,  whereas  the  matantaneous  descriptions  does. 
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